The Asakura-Oosawa model in the protein hmit: the role of many-body interactions 
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We study the Asakura-Oosawa model in the "protein limit", where the penetrable sphere radius 
Rao is much greater than the hard sphere radius Rc- The phase behaviour and structure calcu- 
lated with a full-many body treatment show important qualitative differences when compared to 
a description based on pair potentials alone. The overall effect of the many-body interactions is 
repulsive. 



I. INTRODUCTION 

The Asakura Oosawa (AO) model0], also known as the 
penetrable hard sphere modelQ, was first introduced al- 
most 50 years ago to describe depletion effects in colloid- 
polymer mixtures. The colloids are modelled as hard 
spheres (HS) of radius Rc, and the "polymers", or AO 
particles, as penetrable hard spheres (PHS), which inter- 
act as ideal particles with each other, but as HS of radius 
Rao with respect to the colloids. In spite of its simplic- 
ity, this model has been instrumental in understanding 
the phase-behaviour of polymer-colloid mixtures. For ex- 
ample, increasing the polymer concentration can lead to 
a fluid-fluid or fluid-solid demixing of the colloidal parti- 
cles. The origin of this phase transition arises from the 
depletion effect, which is easily illustrated by calculating 
the free-energy of two HS in a bath of PHS AO parti- 
cles. Each HS excludes a volume |7r(i?c -f Rao)^ from 
the PHS particles, but when two HS approach, their ex- 
clusion volumes overlap, resulting in more free volume 
available for the PHS spheres. This translates into an ef- 
fective depletion pair potential between the two particles 
of the form|l|3 



in the range (7c < r < 2a cp, where (Jcp 
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and ZAO is the fugacity of the AO PHS particles. The 
AO potential is always attractive, with a well-depth that 
increases with the fugacity (or number density) of AO 
particles. Simple geometrical arguments show that for 
size ratios q = Rao/Rc < 0.1547 there are no higher 
order many-body interactions beyond the effective pair 
potential. For larger q, however, multiple overlap of the 
depletion zones can occur, leading to many-body inter- 
actions. As with many soft matter systems, these are 
often difficult to calculate, and so a common approxima- 
tion is to ignore them, and treat a system in the effective 
pair-potential approximation 3, ^] . For the AO model, 
this works surprisingly well for descriptions of phase be- 
haviour up to size-ratios of q 
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In this paper we study the AO model for q = 
Rao/Rc » 1- For polymer-colloid systems this is of- 
ten called the nano-particle or protein limit, because 
small particles such as proteins are needed to achieve 
the large size-ratios. Clearly, a pair potential picture 
should break down for large enough g, where many-body 
effects are expected to dominate. For ideal polymers in 
the limit g < 1, the effect of many-body interaction in 
polymer and colloid mixture has been studied by Meijer 
and Frenkel'sl, who found that these interactions sta- 
bilise the liquid phase. Recent work has shown that for 
q >> 1, many-body effects can qualitatively affect the 
phase behaviour[3, H, ■ 



The AO model was originally developed for size-ratios 
q < 1, where Rao taken to be the radius of gyration Rg of 
ideal polymers 1, 2J. For g >> 1 this simple mapping no 
longer holds, although the AO model can still be mapped 
onto a model for idealpolymers by correctly defining an 
effective Rao radiusQ- However, our goal here is not 
so much to study colloid-polymer mixtures, but rather 
to investigate the effect of many-body interactions on a 
well-defined system. The AO model has the particular 
advantage that the effective pair interaction Vao(^) is 
exactly known. By directly calculating the phase be- 
haviour and structure of a two-component AO model, 
and comparing it to an effective one-component model 
with the AO pair potential, we can systematically study 
the effect of many-body interactions. The insight gained 
from this well characterised system should increase our 
appreciation of the complexity of many-body effects in 
soft matter systems. 



The paper is organised as follows: In section 2 we cal- 
culate the phase-behaviour of the AO model, using Monte 
Carlo (MC) simulations and several simple theories. In 
section 3 we describe the pair correlation functions, and 
also the effective colloid-colloid structure factors. Finally, 
we discuss our results in section 4. 
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FIG. 1: Fluid-fluid binodals for size-ratios q = 3,5,8. Plot- 
ted are the absolute packing fractions r/AO ~ ^t^PaoRao s^nd 
r]c = ^npcRl- Also shown are the binodals calculated with 
the free volume theory of ref.|l3l| and their corresponding crit- 
ical points (asterisks). Note that the critical colloid packing 
fraction tends to zero with increasing q. The free volume and 
virial theory binodals for q — 8 can not be distinguished since 
both theories converge in the limit of large q. 

II. PHASE BEHAVIOUR 
A. Monte Carlo Simulations 

We performed Gibbs ensemble Monte Carlo [loj sim- 
ulations in the semi-grand ensemble, where colloids are 
treated canonically and the PHS in the grand canoni- 
cal ensemble, for three size ratios q — 3,5, 8. The total 
number of colloidal particles {N = 108) is hence fixed, al- 
though they can exchange boxes. The chemical potential 
or fugacity of the PHS is kept constant by the usual grand 
canonical MC insertion and extraction moves Fur- 
ther details of the method can be found in refs. [ill IT^ . 
This setup is equivalent to an AO mixture in osmotic 
equilibrium with a reservoir of only PHS particles[l^. 

Results for the binodals are shown in Fig. For 
increasing q, the critical colloid packing fraction 77c = 
^npcR^ tends to zero, while the PHS packing fraction 
Vao = ^''^PAoR-AO increases. It is also instructive to 
compare the binodals on a log-log plot, shown in Fig.|21 
which emphasises parts of the binodals further from the 
critical point. For example, we see that the binodals cross 
at very low rjc- 

It is often convenient to plot the phase diagrams in the 
semi-grand ensemble, where the colloids are treated in 
the canonical ensemble and the PHS in the grand canon- 
ical ensemble. This is equivalent to setting up an AO 
mixture in osmotic equilibrium with a reservoir of only 
PHS particles p^. Because the PHS particles are ideal, 
the fugacity zao = exp(/3/iyio) = Pao- results in 

this representation are shown in Fig. |3 

The binodals in Figs. El are quite different to those 



FIG. 2: Same as in Fig.Q but on a log-log scale. This empha- 
sizes a different part of the binodals, showing that the relative 
agreement of the free-volume theory with the simulations is 
about the same for all size-ratios. Again, the critical points 
of the free volume binodals are shown as asterisks. Here we 
also include results for q = 1.05, taken from ref. [l^. 
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FIG. 3: Fluid-fluid binodals for size-ratios q = 3,5, 8. The x- 
axis denotes the colloid packing fraction 77c and the y-axis the 
packing fraction rjAo = ^''^PaqRao of a reservoir of pure PHS 
particles at the same chemical potential as the two-component 
AO system. Also shown are binodals calculated with the free 
volume theory of ref. [l^ . which agree very well with the sim- 
ulations. 



of the AO model for g < 1, shown, for example, in 
refs. 6, 13]. First of all, they are at much lower packing 
fractions of the colloids. In fact, the critical colloid pack- 
ing fraction tends to zero for increasing q. Moreover, 
in the reservoir representation, the binodals are much 
narrower, an effect that becomes more pronounced with 
increasing q. Both effects are similar to what has been 
found for mixtures of ideal polymers and colloids 

The physical origin of the low critical colloid packing 
fraction is the large cross-interaction between the two dif- 



ferent species. Perhaps a simpler way of viewing this is to 
first imagine a binary HS mixture of particles of equal di- 
ameter (Tc, but with a cross diameter ai2 > Uc- As ai2 in- 
creases, the system will phase separate at lower and lower 
volume fractions of the HS particles. A phase-diagram 
like that shown in Fig.^would be symmetric in the pack- 
ing fractions of the two species. However, one could also 
choose to represent the packing fraction of one of the two 
species, say species 1, as 771 = ^'Kpiai2- For large enough 
q the phase diagram would then closely resemble that of 
Fig. n since the packing fraction of species 1 would be 
so low that it could be replaced by an AO PHS parti- 
cle without significantly affecting the phase behaviour. 
This way of obtaining the large q limit of the AO model 
also suggests that a simple virial theory should become 
increasingly accurate as q increases and the critical pack- 
ing fractions decrease. The next section will show that 
such a virial theory works very well indeed [l^. 




FIG. 4: The limiting behaviour of the critical point, given by 
Eqs. @ and 0, is compared to the full free- volume theory 
calculations. The differences rapidly decrease with increasing 



B. Two-component free volume and virial theories 



A general Helmholtz free-energy for the two- 
component AO model can be written as: 



F{N,,Nao,V) 
V 



f = f"HPc) + r:^0iPp)+fc~A0iPc,PA0) 

(2) 

where the colloids are treated as hard-spheres, and the 
AO PHS as ideal particles. We suppress the dependence 
on temperature T, since our model is athermal. 

Whereas f ^^ ( Pr ) and Iao^Pao) are well 
understood|l7l|. less is known about the fc-Ao{Pc,PAo) 
term. In the free-volume theory of Lekkerkerker et 
«^-E3> fc-AO is approximated as 



fc-AO = PAO^iVc,q)- 



(3) 



In other words, terms proportional to and higher 
in fc-AO are ignored. uj{rjc,q) can then be interpreted 
as the free energy of inserting a single AO PHS particle 
into a bath of HS colloids at packing fraction rjc- In their 
classic paper, Lekkerkerker et aL[l3j calculated uj{ric,q) 
from scaled particle theory: 



^iVc, q) = - ln(l - rjc) + A 



Vc 



l+Vc 



(4) 
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l+Vc 



c 



l+Vc 



In 



where A ^ 3q + 3q^ + q^ , B = ^q^ + 3q^ , and C ^ 3q- 
the nomenclature of ref. uj{r]c,q) = —\n{a{r]c,q)) 

with a{ric,q) the so-called free-volume fraction, which 
defines the particle packing fraction r]^Q in a reser- 
voir at the same chemical potential as the mixture i.e. 
TjAo — ctiVci <i)Vao- ^s q increases, so does the rela- 
tive strength of the cross interaction, leading to lower 
and lower critical colloid packing fractions. It therefore 



makes sense to expand //TXn in powers of 



rfree 
Jc-AO 



PAO {Vc{l + qf + 0{ril)) 



(5) 



The leading term is proportional to the second cross- 
virial coefficient. In other words, for low rjc (and for 
any pAo), free- volume theory reduces, as expected, to a 
simple virial theory. 

For the virial theory defined by Eq. ||2Jl and the leading 
term of Eq. jSJ, the large q limit of the critical points can 
be derived: 



lim rf^ 



hm 77$£* 



1 



lim 



Vao 



{l + qf 
q' 

a+q)' ' 

exp(l)5^ 
(1 + 9)3 



1 

1 

- exp(l). 



(6) 
(7) 
(8) 



Since r/^*"** ^ as g ^ 00, it follows that free- volume the- 
ory shows the same limiting behaviour, which is demon- 
strated in Fig. ^ 

A comparison with Figs.^-ISlshows that, even though 
higher order pAo effects in fc-AO are ignored, free- 
volume theory agrees remarkably well with simulations, 
just as was found earlier for q « 10, E|. The main 
deviations are found near the critical point; these can 
partially be ascribed to the fact that free-volume theory 
is a mean-field theory, with the wrong critical exponents 
etc., which normally leads to more rounded binodals. 
Nevertheless, the simulations are quite consistent with 
the limiting behaviour for the critical points derived in 
Eqs. ® - ©. 

In Fig. n we also show the binodal for q ~ 8 obtained 
from the simple cross- virial theory. The differences with 
free-volume theory are not visible on the scale of the 
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graph. For smaller q the virial theory is not quite as 
good as the free-volume theory, but it still provides a 
semi-quantitative description of the binodals, suggesting 
that the basic physics can be understood at this simpler 
level. 

The good agreement between simulations and free- 
volume theory for the fluid-fluid binodals suggests that 
we can use the latter to estimate the position of the triple 
point. This occurs when the gas, liquid and solid are in 
equilibrium, and an easy way to estimate the location 
of the triple point is to set up an equilibrium between 
the gas phase branch of the binodal, and a HS fluid at 
the freezing transition. This results in an approximate 
position for the triple point at r^c — 0.494 and: 

^^o"'''«^/3^coe.a,V«6.12q3, (9) 

where the reduced coexistence pressure of a HS fluid at 
freezing, /3Pcoea;Crj?, is known from simulations 10]. In the 
limit of large q it is virtually impossible to fit any PHS 
spheres into the colloidal crystal, while the gas-phase bin- 
odals are at extremely low r^c! we therefore expect this 
relationship to become asymptotically exact for large size 
ratios. In fact, Eq. © gives a good prediction for the 
free-volume triple point for all size ratios where a triple 
point exist, i.e. even for q < 1. For example, at g = 0.8 
free volume theory gives yy^Q = 3.13 while Eq. © would 
predict rfj^Q — 3.17. This analysis shows that, in the 
semi-grand ensemble, the triple point moves to extremely 
large values of ?7^o, compared to the critical point. 

C. One-component theory with pair potentials 

It is instructive to compare the results for the two- 
component AO model with those of an effective one- 
component model. As mentioned in the introduction, 
the AO PHS particles can be integrated out to derive an 
exact pair potential, valid for all g, and given by Eq. JQ. 
For q < 0.1547 this leads to an exact description of the 
system, but for larger q, many-body interactions must 
be invoked. Nevertheless, for practical reasons, the pair 
approximation is often used in soft matter physics. In 
many situations this works well, but here we expect it to 
break down as q increases. 

For large q, the pair potential of Eq. ^ becomes very 
long-ranged with respect to the colloidal diameter CTc- 
Therefore mean-field theory, for which the free-energy 
takes the form 

y = fc [Pc) + ^Pc J drVAoir). (10) 

should become asymptotically exact. For this simple 
"van der Waals limit" , the critical colloid packing fraction 
is always given by rj^"* ss 0.13, independent of potential 
details. Here we work in the semi-grand ensemble, where 
the effective pair potential picture has a consistent statis- 
tical mechanical interpretation0,0|. In Fig. [S] the bin- 



o 



o 

ti - 





q =3 Sim 
£7=3 MF 


1 - 

0- 




A £7 =5 sim 
£7=5 MF 

o n —ft c/Vn 
^ — o oil 1 1 


-1 - 


V 


£7=8 MF 


-2- 


, 




-3- 






/I 







0.0 0.1 0.2 0.3 0.4 



% 

FIG. 5: A comparison of the binodals for the two-component 
AO model, calculated by GEMC simulations, with an effec- 
tive one-component picture, calculated with the pair potential 
VAo{r) and Eq. IIOII . Note the qualitative differences in the 
limit of large q: For the one-component picture ryj;'^'' tends to 
0.13, and rfj^Q^^ tends to 0, whereas the full two-component 
picture shows the opposite behaviour: Ty^"^** tends to zero, and 
'7ao^'' tends to a constant. These trends are depicted by the 
arrows. The differences are due to the effects of many-body 
interactions. 



odals from Eq. H10|l are compared to the full 2-component 
simulations. 

Besides the obvious quantitative differences for the lo- 
cation of the critical points (note that the y-axis is log- 
arithmic!), there are important qualitative differences as 
well. For example, the two-component binodals are much 
more narrow. Furthermore, the critical points show op- 
posite scaling behaviour with increasing size ratio q: For 
the two-component model rfj^Q^*' , tends to a constant and 
rjc tends to zero, whereas for the one-component model 
^AO** fsnds to zero and r]c tends to a constant. 

The dominant effect of the many-body interactions ap- 
pears to be repulsive, since phase-separation occurs at a 
much higher packing fraction of the AO PHS particles. 
In fact, the three-body interactions have already been 
calculated for the AO model ||2j|; they are repulsive for 
all geometries. This would seem consistent with the over- 
all effect of the many-body interactions. However, this 
interpretation is most likely too naive. For example, we 
expect that the fourth order term is attractive again, and 
that the series oscillates, as was recently found in simu- 
lations of a self avoiding walk polymer system[22. In 
general, the sign of a many-body interaction can vary 
in a complex way with coordinates. A good example is 
given by the 3-body HS depletion interactions calculated 
in ref. 20]. Furthermore, it has recently been shown that 
the effects of many-body interactions in a mixture of in- 
teracting polymers and colloids are attractiveH,|22, the 
opposite of what we find for the AO model. These ex- 
amples suggest that it is generally quite difhcult to make 
simple predictions regarding the effect of many-body in- 
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teractions on phase-behaviour. 



III. PAIR STRUCTURE 
A. Radial distribution functions 

Given the unusual phase-behaviour of the AO model in 
the protein limit, it should be interesting investigate the 
effect of many-body interactions on the pair structure. 
To that end we performed Monte Carlo simulations of 
the colloid-colloid and colloid- AO pair correlations. Rep- 
resentative examples are shown in Fig. |B1 Both gcci^) 
and gc-Aoif) have fairly weak structure, with only one 
main peak. The reasons for this are two fold: 1) The 
colloid packing fraction is very low. 2) More generally, 
long ranged interactions lead to less sharply peaked radial 
distribution functions We also compared the simula- 
tions to pair correlations derived from the test-particle 
route for a recently developed fundamental measure the- 
ory (FMT) density functional theory (DFT)0i]. This 
DFT shows the same phase-behaviour as the free volume 
theory, which we have shown to be very accurate for the 
AO model. We therefore expect the DFT to be rather 
good in this limit, although correlations are often a more 
sensitive probe of a DFT than phase behaviour is. Re- 
sults are compared in Fig.^jto the simulations. For lower 
r]AO the agreement is very good, but for higher rjAo some 
deviations are found. This may partially be because the 
state point for the highest 77^10, {Vc — 0.004, 77^0 = 0.6) 
is very close to the free- volume critical point (which lies 
at ijc — 0.0046, TjAo — 0.586), while the critical point for 
the simulations is further away. 

On might expect that, due to the low colloid densi- 
ties, the correlation functions should be rather straight- 
forward to calculate with integral equation methods[l7l|. 
In Fig.[7|we compare, for one state-point, some represen- 
tative results from the Percus Yevick (PY) and Hypernet- 
ted Chain (HNC) approximations [T^ . PY systematically 
underestimates the peaks, a general effect that becomes 
more pronounced with increasing rjAo : f^nd which is sim- 
ilar to what happens for binary HS mixtures^- HNC 
appears to be more accurate, but it suffers from a non- 
solution line, a mathematical artifact where no solutions 
are found. This non-solution line occurs well before the 
expected spinodal, and makes HNC less useful. Finally 
we note that direct functional differentiation of the FMT 
DFT in the bulk leads to PY correlation functions. The 
test-particle route we apply here is thought to generally 
be more reliable. 



B. Structure factors 

For scattering experiments, a more useful measure of 
the pair correlations is given by the colloid-colloid struc- 
ture factor S{k). These are shown in Fig.|Hlfor some of 
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FIG. 6: The colloid-colloid and coUoid-PHS radial distribu- 
tion functions for g = 5 at several state-points. We compare 
direct simulations with results from the test-particle route of 
a fundamental measure theory DFt|23|. For lower packing 
fractions the agreement is good, but when the DFT binodal 
is approached, differences become larger. 




FIG. 7: A comparison between DFT, HNC, and PY ap- 
proaches to the colloid-colloid structure, and results from sim- 
ulations. 
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FIG. 8: Structure factors from simulations for some of the 
same state points as Fig. |^ Inset: the first isosbestic point 
where the structure factor is invariant to changes in rjAO- 



FIG. 9: Structure factors from HNC for different values of 
q. The state points are chosen to give a similar value of 5(0). 
The main effect of increasing q is to decrease the value of fcuc 
where S{k) begins to increase towards its maximum at S(0). 



the same state points as in Fig. El The S{k) shows vir- 
tually no structure, except for a maximum at S{Q). This 
behaviour seems to be generic for the large q limit of the 
AO model. Another interesting observation is the occur- 
rence of an isosbestic point kisos, a value of k where S{k) 
is invariant to changes in r]Ao- A recent theorvp3| for 
isosbestic points predicts that kisosO'c decreases with in- 
creasing range '25'. The theory predicts that kisosCc ~ tt 
for very short ranged potentials, and that for the AO 
model it can be approximated as: 



hsosOc « V (1 + 0.42(7/2 + (0.42g)Vl2) (11) 



Even though this theory was derived for the small q limit 
of a one-component model, it still appears to be semi- 
quantitative here, since we find an isosbestic point at 
kisos = 1.3±0.05 and Eq. ifTT)) would predict kisos = 1-30. 
For other q's we also found good agreement. This fur- 
ther suggests that one experimental signature of the long- 
ranged nature of the AO model is the late upturn of S{k) 
toward a maximum at 5(0). This is illustrated in Fig. 13 
We note that in a recent paper, Tuinier and Brulet(26i| 
have shown similar behaviour from a one-component cal- 
culation with only effective pair potentials. They also 
performed small angle neutron scattering experiments on 
lysozyme-polysaccharide mixtures, measuring structure 
factors S{k) that appear quite similar in shape to those 
predicted here for the AO model. Since these authors 
obtained similar qualitative behaviour to our simulations 
from a one-component theory, the qualitative behaviour 
of S{k) that we observe is most likely caused by the long- 
ranged nature of the pair potentials, and not so much by 
the many-body character of the interactions. 



C. Effective colloid-colloid pair potentials 

The extremely low packing fractions of rjc studied here 
might suggest that one could approximate the pair cor- 
relation function by its zero density limit: 

lim gUr) = exp [-I3Vao (r)] . (12) 

However, this is not correct, as suggested by the large 
differences between the binodals calculated only with 
yAoi^) and those calculated for the full two-component 
AO model. We demonstrate this explicitly in Fig. 1101 
where gcc{f) at the same ry^Q (which means the same 
Vko (?")), but for different rj^ was calculated with the 
FMT DFT approach. At -q^ = 0, eq. lO is of course 
exactly obeyed, but this no longer holds for the other 
values of rjc- The effect of the many-body interactions 
lowers the peak of g{r), which is consistent with the be- 
haviour of the binodals, where the effect of the many- 
body interactions is to reduce the cohesion between the 
HS colloids. 

A general theorem states that for any g{r) and density 
p, there exists a unique pair potential Wg(r; p) which will 
reproduce that (7(r), regardless of the underlying many- 
body interactions [23. We inverted the gcdr) at two 
state-points with HNC and PY inversions. The results 
are shown in Fig. 1111 The two inversion methods give 
very similar results (especially at low polymer packing 
fractions, where both inversions can not be distinguished 
on the scale of the graph), suggesting that the potential 
obtained is indeed close to the true Vg{r;pc). Note that 
Vg{r; pc) is different from the simple potential of mean 
force —ln{gcc{r)), which shows that correlation effects 
are important here. As expected, Vg{r; pc) is less attrac- 
tive than the bare pair potential Vao (t) ■ It is also slightly 
longer ranged. If there were no many-body forces, then 
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FIG. 10: Radial distribution functions for systems with differ- 
ent rjc- For eadi ric the reservoir packing fraction was rj^o ~ 
0.964792, which corresponds to rjAO = 0.4 at rjc = 0.00402. 



Vg{r; Pc) would be equal to Vao(^) at all densities. The 
differences can therefore be attributed to many-body in- 
teractions, whose overall effect is to weaken the effective 
pair potential. 

At each density one could use Vg{r;p) to extract the 
colloid-colloid pair correlations, and also the osmotic 
pressure through the compressibility route. However, in 
general a large density dependence of the effective pair 
potential also implies difficulties like the representability 
issues discussed in ref . , making a pairwise description 
less useful. 



D. Effective PHS-PHS pair potentials 

In the limit of large q it might seem tempting to try 
the opposite of the usual Asakura-Oosawa strategy, and 
integrate out the HS colloids instead of the PHS parti- 
cles. This can easily be done, and the potential between 
two isolated PHS particles takes the same functional form 
as in Eq. QJ, but with colloid and PHS AO parameters 
switched. One difference is that Vao('') is now relevant 
down to r = 0, because there is no HS repulsion to re- 
strict it to r > (Tc: as is the case for the colloids. However, 
a description based only on the effective pair potential 
immediately leads to problems. The lack of hard-core re- 
pulsion means that the effective one-component system 
of PHS particles with an AO type pairwise attraction 
falls into the class of catastrophic potentials defined by 
Ruelle28], for which there is no thermodynamic limit. 
This does not mean that one cannot derive a consis- 
tent thermodynamics by integrating out the HS colloids. 
Rather, because the PHS particles can overlap so easily, 
many-body effects are always very important, especially 
at small r, and are necessary to stabilise the effective PHS 
system. A similar situation was recently found when in- 
tegrating out internal degrees of freedom for a solution of 
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FIG. 11: Inversions of the colloid-colloid pair correlation 
functions for two state-points: up: 77c = 0.0402, t/ao = 0.1 
and down: rjc — 0.0402, rjAO = 0.6. For each we show the re- 
sults of a PY and an HNC inversion, as well as the potential 
of mean force, — \n{gcc{r)) and the bare two-body potential 
VAo{r)- HNC and PY are indistinguishable for the lower 
polymer packing fraction (upper figure). 



polymers in a p oor solvent 29], as well as for a Gaussian 
Core model Both examples can lead to catastrophic 
pair potentials, even though the underlying many-body 
system is stable. 



IV. CONCLUSIONS 

In conclusion then, we have shown that in the so-called 
protein limit, where g >> 1, the behaviour of the two- 
component AO model differs significantly from a descrip- 
tion based on an effective pair potential description alone. 
We summarise our main results: 

• In the limit of large q, the AO binodals move to 
lower and lower colloid packing fractions r]c, and at 
the critical point, the AO particle packing fraction 
rjAO tends to a constant. This is opposite to what is 
found from a pair potential description alone, and 
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so this phase behaviour can be ascribed to many- 
body interactions. 

• The free- volume theory of Lekkerkerker et ol.fl^ 
works remarkably well for the phase-behaviour. It 
reduces to a simple virial theory in the large q limit, 
allowing us to extract limiting values for the critical 
points which are consistent with the simulations. 

• The colloid-colloid pair-correlation functions show 
very weak structure. Even at very low packing frac- 
tions rjc, gcc{f) is not well described by its zero- 
density limit exp[—pVAo{f )]- The pair correlations 
are well described by a recent fundamental measure 
theory density functional theory ,2^], at least if one 
is not close to the critical point. 

• The overall effect of the many-body interactions is 
repulsive, as seen in the phase behaviour, and also 
in the structure. 

• There are clear signatures of the long-ranged inter- 
actions in the structure factors S{k). In particular, 
the value of kac where S{k) begins to rise to its 
maximum at fc = 0, decreases with increasing q. 

• A description based on effective pair potentials be- 
tween the PHS particles, derived by integrating out 
the smaller HS colloids, leads to catastrophic sys- 
tems with no thermodynamic limit. 



Our aim in this paper was to study the effects of many- 
body interactions in a well-defined model system. Al- 
though it would be tempting to extract some more gen- 
eral insights about the role of many-body interactions 
in soft matter systems, this is not so easy to do. On 
the one hand, we can make predictions about the be- 
haviour of a related many-body system, namely a mix- 
ture of ideal polymers and HS colloids in the limit of 
large g0, where we expect some similar trends. But 
on the other hand, if the ideal polymers are replaced 
by interacting ones, the behaviour changes: for example, 
the critical colloid packing fraction is almost constant, 
and the overall effect of the many-body interactions is 
attractive instead of repulsive^. It is clearly not always 
easy to predict the effect of the many-body interactions 
a-priori. We conclude from this that coarse-graining a 
soft-matter system to a representation where many-body 
interactions are important when compared to the pair in- 
teraction may not always be a very fruitful way forward. 
Sometimes it may be easier to treat the original system 
without this coarse-graining step. 
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